Thermodynamics of the vortex liquid in heavy ion-irradiated superconductors 
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It is shown that the large effect of heavy ion-irradiation on the thermodynamical properties of 
the anisotropic superconductor YBa2Cu307_i extends well into the superconducting fluctuation 
regime. The presence of the induced amorphous columnar defects shifts the specific heat maximum 
at the normal-to-superconducting transition. This effect is similar to that recently put into evidence 
in cubic K^Bai-^BiOa [x — 0.35). In both compounds, vortex pinning manifests itself as a sharp 
angular dependence of the equilibrium torque. In YBa2Cu307_<5, pinning by the defects appears 
at the temperature T™ ax of the specific heat maximum, well above the magnetic irreversibility line 
Ti r r(H). In isotropic K^Bai-ajBiOs, the onset of the pinning-related torque anomaly tracks the 
onset of the specific heat anomaly and the irreversibility line. In YBa2Cu307_,5, fluctuations of the 
amplitude of the order parameter (and not vortex line wandering) are ultimately responsible for the 
vanishing of pinning. In K^Bai-^BiOs, vortex pinning disappears only at the superconducting-to- 
normal transition. The results indicate that in both compounds, the pinning energy at the "Bose 
glass" transition is large with respect to the total free energy gain in the superconducting state. By 
implication, the mechanism of this latter transition should be reconsidered. 

PACS numbers: 74.25.Bt, 74.25. Op,74.25.Qt,74.70.+k 



I. INTRODUCTION 

When heated above their irreversibility line, disor- 
dered type II superconductors undergo a transition from 
a "truly superconducting" ensemble of localized vortex 
lines to a "vortex liquid" of diffusing lines. This transi- 
tion is most commonly described in terms of a thermal 
"depinning" of vortex lines from material defects. Typ- 
ically, thermal wandering of the vortex lines from the 
defects is thought to become increasingly important as 
the temperature is raised, until, at the transition, the 
free energy gain obtained from vortex localization on the 
defects has dropped to ~ fcsT. 

The case where the defects are columnar amorphous 
tracks introduced by swift heavy-ion irradiation has at- 
tracted much attention, not in the least because the prob- 
lem becomes particularly tractable theoretically. Using 
the formal analogy between flux lines and 2D bosons 
in a static disorder potential, Nelson and Vinokuni cal- 
culated the magnitude of thermal positional excursions, 
concomitant pinning energies, and the resultant phase 
diagram. This consists of a low temperature disordered 
"Bose glass" of localized vortices that gives way to the 



vortex liquid at the temperature Tbg(B) [or induction 
Bbg{T)\. Model descriptions of the Bose glass-to- liquid 
transition nearly exclusively rely on the total pinning en- 
ergy near Tbg{B) being small. This can be either be- 
cause line wandering reduces the single-vortex pinning 
encrgy^^ or because vortices vastly outnumber the de- 
fects. In both cases, the pinning energy gain is only a 
small perturbation to the total free energy of the system, 
and the position of the Bosc-glass transition line can be 
obtained as a shifti*2ii^i& of the (first order) vortex lat- 
tice melting line of the pristine materialise This ap- 
proach has had some success in explaining the observed 
increase of Bbg with defect density rid in YBa2Cu307_5, 



1/2 1/2 
n d 



(with 



where it was found that 8B B g /9T ~ 1+A$ t 
3>o = h/2e the flux quantum)i2i4 

Simultaneously, experiments abound indicating that 
the pinning energy near the Bose-glass transition in 
cuprate superconductors is not small. Reversible mag- 
netization measurements on heavy-ion irradiated layered 
Bi2Sr2CaCu20§ reveal a large contribution of columnar 
defect pinning to the free energypi°*ii An effect of the 
ion tracks was measured up into the fluctuation critical 
regimes*, and interpreted in terms of local T c -variations 
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induced by the defectsii2ii£ Recent measurements on cu- 
bic (K,Ba)BiC>3 have shown that heavy-ion irradiation 
even affects the specific heat C p : the temperature at 
which the C v jump occurs, signaling the transition to 
the superconducting state, was found to shift upward 
with increasing defect density rid, and depends on the 
angle between magnetic field and the track direction^&i^ 
When the magnetic field is aligned with the tracks, su- 
perconductivity is enhanced, when the magnetic field is 
turned away from the tracks, one recovers the behav- 
ior of the pristine crystal. As for heavy-ion irradiated 
YBa^CusOy-i, transport experiments show that in the 
vortex liquid, the resistivity remains exponentially small 
with respect to that of the pristine m aterial^^ mSi 
More strikingly, the experiments of Refs. Ilfill7ll8l reveal 
an angular dependence of the resistivity, related to vor- 
tex pinning by the tracks, that persists up to resistance 
levels that are 90 percent of that in the normal state. In 
other words, columnar defects affect vortices in the vor- 
tex liquid, and up into the fluctuation paraconductivity 
regime. However, there are no thermodynamic measure- 
ments assessing the importance of the pinning energy in 
this material. 

In order to establish the magnitude of the contri- 
bution of vortex pinning by columnar defects to the 
free energy, we have performed measurements of the 
specific heat and the reversible torque on heavy ion- 
irradiated single crystalline YBa 2 Cu307_(5 fSections lII CI 
and 111 Bfl . The data are compared to previous specific 
heat resultsiiii^ and new torque measurements on single 
crystalline K^Bai-^BiOs. The main difference between 
the two materials lies in their Ginzburg number Gi = 
i [fc B T c /27T£ £o (0)£(0)] 2 . Here e (T) = $ 2 /47r Mo A 2 b (T) 
is the vortex energy scale, X a b{T) is the penetration 
depth for currents running perpendicularly to the crys- 
talline anisotropy axis, and £(T) the coherence length (for 
YBa 2 Cu307_«5, £ = £ &, the a6-plane coherence length). 
Since Gi for optimally doped YBa2Cu30y_5 is two or- 
ders of magnitude larger than that of K^Bai-^BiOs (see 
Tabled, it is conceivable that thermal fluctuations wipe 
out any strong effect of columnar defects in the vortex 
liquid phase. On the contrary, we find that the reduced 
temperature T™® X {H)/T C at which the specific heat is 
maximum in YBa2Cu307_,5 unambiguously shifts up- 
ward with increasing columnar defect density rid, as it 
does in K^Bai-^BiOs. This upward shift can only be 
accounted for by a large (i.e. not perturbatively small) 
contribution of pinning by the columnar defects to the 
free energy. This is in contradiction to the assumptions 
commonly made in estimating the Bose-glass transition 
HneiSi&aiSiS In Section UlI Bl we estimate the mean-field 
pinning energy contribution required for the specific-heat 
shift. 

We also obtain the experimental pinning energy di- 
rectly from reversible torque measurements (Section 
IIII All . It turns out that the field- and temperature de- 
pendence of the pinning energy, rather surprisingly, scales 
with the parameter Q = (1- &)(l-t 2 ) 1 /3(^)-2/3 G - i -i/3 ) 



suggesting that a development of the Ginzburg-Landau 
free energy functional in terms of Lowest Landau level 
(LLL) eigenfunctions is an appropriate starting point for 
a model description. 19 i 20 i 21 i 22 i 23 i 24 

Here, t = T/T* IF with T C MF the mean-field transi- 
tion temperature, and b = B/B C 2(T) with B C 2(T) = 
<&o/27r£ 2 (T) the upper critical field. An assessment of the 
experimental result shows that in YBa2Cu307_a, fluctu- 
ations of the order parameter amplitude lower the pinning 
energy with respect to the expected mean-field value. In 
K K Bai_2;Bi03, with small Gi ~ 1CP 5 , thermal fluctu- 
ations are unimportant and pinning subsists up to the 
superconducting-to-normal state transition. 

II. EXPERIMENTS AND RESULTS 

A. Samples 

Experiments were done on a series of untwinncd and 
lightly twinned YBa2Cu307_a single crystals, grown 
by the flux method in Au crucibles, and subsequently 
annealed in oxygen in Pt tubes^ Measurements on 
K^Bai-^BiOs were made on a crystal with x = 0.35, 
grown by electrocrystallization. The crystals were irra- 
diated during different runs at the Grand Accelerateur 
National d'lons Lourds (GANIL) in Caen, France. The 
YBa2Cu307_,5 crystals were irradiated with 5.8 GeV Pb 
ions; for all but one crystal, the beam was aligned parallel 
to the c-axis. The final crystal (hereafter referred to as 
"Y30" ) was at an angle of 30° with respect to the c axis, 
to a fluence of 1 x 10 11 ions cm~ 2 , which corresponds to 
a matching field = <&$nd = 2 T. The K .35Bao.65Bi03 
crystal was irradiated with 7.2 GeV Ta ions. The irradi- 
ation produced continuous amorphous columnar defects 
of radius cq « 3.5 nm. The irreversibility line Ti rr (H) 
of all crystals was measured as the onset temperature 
of the third harmonic of the ac transmittivity Tnsr^^ 
with the DC field aligned parallel to the defect direction. 
The irreversibility line of the "Y30" and Ko.35Bao.6.sBi03 



YBa 2 Cu 3 07-i K .3sBao.65BiO 3 



A A ab (0) = 120 nm A(0) = 220 nm 

£ Cab(0) = 1.4 nm £(0) = 3.8 nm 

eo(0) 1.7 x 10" 11 Jm^ 1 5.2 x 10" 11 Jm" 1 

Gi 2 x 10~ 3 1 x 10~ 5 

S = \ab/K = tic/tab 0.14 1 



TABLE I: Superconducting parameters for the studied com- 
pounds. 
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FIG. 1: YBa 2 Cu 3 07_,5 : Torque signal for rotation angles 
9 close to the irradiation direction 8d, for different angles ■& 
between the plane of rotation of the magnetic field and the 
defects (H = 20 kOe, T = 88 K). 



crystals was also obtained from torque magnctomctry 
(see Table |nj. 
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FIG. 2: YBa 2 Cu 3 7 -i : Slope dV/dO of the torque signal 
for field along the c axis, divided by field. Lines are linear 
fits extrapolating to T C (H). Inset (b): line energy £0(2^) ob- 
tained by dividing out the field dependence from the data 
from the data in the main panel. Inset (c): Scaled magne- 
tization H~ 1 dT /d9/(tb) 2 ' 3 as function of the LLL parameter 
Q = (I-fe)(l-f 2 ) 1/3 /(rt>) 2/3 Gi 1/3 . 



B. Torque measurements 

Apart from the characterization of the pinning energy, 
torque measurements were also used to obtain the su- 
perconducting parameters of the two compounds under 
study (see Tabled . The measurements on YBa2Cu307_ ( s 
were performed on the same twinned single crystal (Y30) 
as that used in Ref . |2£| The microtorque setup^ was im- 
proved by the adjunction of a secondary 
magnetic field perpendicular to the main one, so that 
the plane in which the field is rotated could be chosen 
arbitrarily. The angle -d between this plane and the de- 
fect direction could be set with a resolution better than 



0.1°. We set •& to zero ±0.1° by maximizing the pinning 
measured by the torque irreversibility along the defect 
direction (Fig. [TJ. 

The mean field upper critical field line B C 2(T) was lo- 
cated using the slope of the equilibrium torque per unit 
volume. For small angles 9 between the applied field and 
the c axis direction, this is given by 



dT/d0 ~ HM±(H), 



(1) 



where Mj_ is the magnetization for the field applied 
along the c-axisi 29 i 30 i 31 A plot of H~ 1 dT/dQ, shown in 
Fig. |2 represents M±(H). For the lower temperatures, 
H~~ 1 dT/d9 depends linearly on T, in agreement with the 
(mean-field) 
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pristine 


= 1 T 


= 2 T 




B^ = 2 T, 30° 




Dimensions (^im 3 ) 


330 x 400 


triangle of base 650, 


430 x 510 


200 x 450 


130 x 337 x 18 


120 x 45 x 30 




x20 


height 530, thickness 20 


x20 


x20 


"Y30" 




Description 


untwinned 


16 TB's || base 


TB's spaced 


twinned 


TB's spaced 


tracks | 




rectangle 




by 10 fim 




by 5 /im 


long dimension 


Tq£ x (H = 0) ( K ) 


93.1 


92.4 


92.1 


91.0 


* 


* 


TL{H = 0) ( K ) 


93.1 


92.5 


92.0 


90.8 


91.3 


31.5 


T C MF ( K ) 


93.1 


93.1 


93.1 


93.1 


92.3 


32.4 



TABLE II: Characteristics of single crystals used in this study. 
All crystals were either untwinned or twinned with a single 
twin boundary (TB) orientation. 
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FIG. 3: YBa2Cu307_^: magnitude of the torque discontinu- 
ity at field alignment with the defect direction 6d- Straight 
lines represent the extrapolation of the torque discontinuity to 
Tk(H). The inset shows the torque curves for both directions 
of rotation at ^lqH = 2 T and T = 0.97T C . It also shows the 
determination of the torque discontinuity To by extrapolation 
of the linear signal to 8 — 0. 



Abrikosov formula 

with (3a = 1.16. According to Eq. @, the zero intercepts 
of H~ l dT/d0 correspond to T C (H), i.e. the B c2 (T) line. 
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FIG. 4: YBa 2 Cu 3 07-6: Comparison of the H c2 {T) line ( □ ), 
TSS*{H) (•), T k {H) (O ), and T lrr (H) (., O ) for the two 
crystals with n d = lx 10 11 tracks cm" 2 {B$ = 2 T). The 
drawn line denotes the locus oi\F n — F s (B)+n t U p (B)\ = 
aksT/V, with V — ea^ab and a — 1 (see text). The 
dotted line shows the same, but choosing V — e^ b and 
a = 0.1. 
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FIG. 5: Ko.ssBao.esBiOs: Amplitude of the torque discon- 
tinuity To for field alignment with the defect direction 9d- 
Straight lines show the extrapolation of the torque step to 
Tk(H). The inset shows the torque curves for both direc- 
tions of rotation and the determination of the torque jump as 
the difference between the linear extrapolations of the torque 
from high angle to 6 — (dotted lines). The applied field 
fi a H = 1.25 T, the temperature T = 0.92T C . 



A plot of the zero intercepts for different fields yields the 
slope dB c2 /dT = -2 TK" 1 and £ ab (0) = 1.37 nm, while 
dividing out the field dependence 1 — B / B c2 (T) gives the 
line energy eo(T). The inset (b) of Fig. |3 shows that 
£q(T) docs not depend on field, as it should. The re- 
sults were checked by plotting the scaled magnetization, 
H~ x dT I dQ I (tb) 2 / 3 as function of the LLL parameter Q, 
see inset (c) of Fig. [21 The magnetization could be scaled 
using the experimental dB c2 /dT = —2 TK _1 and the 
mean-field critical temperature T C MF = 92.3 K. The lin- 
ear portion of the H~ l dT/d6 curves falls in the regime 
Q > 10, well outside the fluctuation-dominated region. 
The S C 2(T)-line, plotted in Fig. [3] the £ a b value deduced 
from it, and de^/dT, which yields the penetration depth 
extrapolated to zero temperature A a b(0) = 120 nm, agree 
well with the results in the literature^ 

As first shown in Ref. l28l and reproduced in the In- 
set of Fig. [21 the equilibrium torque features a step of 
magnitude Tq at the angle at which the field orientation 
coincides with the defect direction 6 = 6d- On both sides 
of the step, the torque signal depends linearly on the field 
angle, allowing Tq to be determined from the extrapola- 
tion of the torque to 0—6d = (Inset to Fig.[3J). The main 
panel of Fig. shows Tq as function of temperature, for 
fields O.5i?0 < fXoH < 1B$, The torque step decreases 
approximately linearly with temperature, with a roughly 
field-independent slope dTo/dT = 65 ± 5 Jm _3 K -1 , and 
vanishes at the onset temperature Tk(H) [or, conversely, 
at the field H k (T)]. The locus of T k (H), as well as the 
irreversibility line Ti rr (H) and the B c2 (T)-line, are dis- 
played in Fig. 0J Additional measurements on samples 
with B = 1 T and B^ = 0.4 T indicate that T k (H) 
depends weakly on the irradiation dose in this range of 
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FIG. 6: Ko.ssBao.esBiOa 

cific anomaly heat onset temperature T<£" set (i/) ( □ ), the 



B,), = 2 T: Comparison of the spe- 

rions 
Cp 



specific heat maximum Tg£ x (H) (•), T k {H) (O ), and 
Ti rr (H) ( O ) . The drawn line denotes the locus of 
\F S (B) + n t U p {B) -F n \ = ak B T/V, with V = afe and 
a = 1 (see text). The dotted line shows the same, but 
choosing V = £ 3 and a = 0.1. 



FIG. 7: YBa2Cu307_5: Normalized specific heat data on a 
pristine crystal (points), and on a crystal with rid = 1 x 10 11 
cm -2 , i.e. B,f, = 2 T (thin lines). In order to compare the 
two data sets, the temperature has been rescaled to t — T/T c , 
and the specific heat to the normal state (phonon) contribu- 
tion at T c . The temperature dependence due to phonons was 
subtracted. The inset shows an enlargement of the specific 
heat curve measured on the pristine crystal in fields of 3 and 
4 T. 



B<f,. 

Concerning Ko.35Bao.6sBi03, a discontinuity in the re- 
versible torque signal, superimposed to some irreversibil- 
ity along the track direction, can be put into evidence 
in a way similar to YBa2Cti307_5. The torque step 
again depends roughly linearly on temperature, with 
slope dTo/dT sa 6 ± 0.5 x 10 2 Jm" 3 !^ 1 (Fig. 0). The 
Tfc(iJ) line is found to lie slightly above Ti rr (H) (Fig- 
ure ini . 

For both compounds, the temperature- and field de- 
pendence of To is found to be parameterized, at all but 
the lowest measuring fields (/io-ff = 1-1 T for Y30, 0.75 
T for K .35Bao. 656103), by the LLL parameter Q [see 
Fig. |S] (a,b)]. All torque step data trace the same curve 
when divided by the free energy density ksTB /$o£(0) 
and plotted vs Q, assuming the mean field transition 
temperature to be T* IF = 93.1 K in YBa 2 Cu 3 07-5, and 
32.4 K in Ko.35Bao.65Bi03. The condition for thermo- 
dynamic quantities to depend on the sole parameter Q 
is that the Landau level splitting ~ 2T^ IF h be greater 
than the strength of fluctuations T™ F (2Gi) 1 / 3 (ht) 2 ^ 3 
[with h = 2tt£ 2 (0) j B/$ = B/T™ F /(dB c2 /dT) T=T Mv = 

(db^ 1 /dt)^li\~ In YBa2Cu307_,5, at the temperatures 
under investigation, this condition is satisfied for fields 
larger than 1 T. The failure of the low-field data to con- 
form to the scaling might thus be due to Landau level 
degeneracy. In Ko.3sBao.65Bi03, however, the LLL con- 
dition is supposed to be satisfied at all fields investi- 
gated. The fact that the low-field data do not conform 
to the LLL-scaling in this compound either, suggest an- 



other origin of the breakdown of scaling. Namely, the 
intrinsic inhomogeneity introduced by the randomly po- 
sitioned columnar defects leads to a spread in local field- 
dependent critical temperatures to which the supercon- 
ductor is sensitive at fields B < B^w&H As a result, the 
effective T c is not the real critical temperature but an 
average quantity (with a lower value) determined by the 
defect distribution. This effect was previously put into 
evidence in heavy- ion irradiated Bi2Sr2CaCu208+<5: at 
fields O.2i?0 < B < Bcf, the magnetization was found to 
approximately follow the LLL scaling relation but with 
an effective critical temperature 2 K lower than the mean- 
field T c of the pristine sample. We suggest that a sim- 
ilar effect is responsible for the downward shift of the 
present 1 T data in Y30, and the 0.75 T data on the 
Ko.35Bao.65BiO-3 crystal, as well as for the difference in 
T C MF deduced for Y30 for field along the tracks and field 
along the c-axis (in the latter case, the field component 
along the defects is smaller). 

C. Specific heat 

Specific heat measurements have been performed on all 
YBa2Cu307_i crystals (see table lTT|) . The same measure- 
ment technique was employed as in Ref . Well-defined 
specific heat anomalies were observed in all crystals, in 
spite of the use of large irradiation doses (up to B^ = 5 
T). 

In the pristine sample, the amplitude of the zero field 
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FIG. 8: Lowest Landau level scaling of the torque (a,b) and 
specific heat data (c-f). (a) Scaled torque jump, measured 
on YBa2Cu307_i crystal Y30, versus LLL parameter Q (see 
text). The parameter T^ IF — 93.1 K. (b) Scaled torque jump, 
measured on the Bao.65Ko.3sBi03 crystal with = 2 T, 
versus Q, with T C MF = 32.4 K. (c) C/T versus Q for irradiated 
YBa 2 Cu 3 07-i with = 2 T. (d) C/T versus Q for pristine 
YBa 2 Cu 3 7 _i. 



superconducting jump was of the order of 4% of the to- 
tal specific heat, attesting to its very high quality. The 
anomaly presents the typical shape of the superconduct- 
ing transition in presence of strong thermal fluctuations 
for H = 0, broadens for increasing H and is shifted to- 
wards lower temperature (see Fig. 0). Sharp vortex lat- 
tice melting peaks are observed for 1 T< H < 6 T (see 
inset of Fig. 

The presence of the amorphous columnar defects re- 
duces the absolute temperature Tq^ x at which the spe- 
cific heat maximum occurs in zero field (see table[p*A 
lowering of the critical temperature after heavy ion irra- 
diation may occur as a result of "self-doping" of the inter- 
column space by O ions expelled from the tracks^^l but 
no such effect was reported for YBa2Cu307_,5- Rather, 
the columns may act by reducing the average zero field 
T c at which long range superconducting order can set 

Turning to the specific heat data in non-zero field, we 
find that the columnar defects increase the (reduced) 
temperature T^ X (H)/T C of the C p maximum (Fig. [J. 
Furthermore, the specific heat curves in magnetic field 
become sharper after irradiation. Apparently, columnar 
defects suppress order parameter fluctuations in a mag- 
netic field, even though the effect in YBa2Cu 3 07_,5 is 
weaker than that previously found in K x Bai_j;Bi03iiii^ 



No vortex lattice melting anomaly is observed in the ir- 
radiated crystals, nor is any other anomaly at, e.g., the 
irreversibility line. 

Figure 0| shows that in YBa2Cu307_,5, the maximum 
of the specific heat systematically coincides with the tem- 
perature Tfc at which the torque signal from vortex pin- 
ning by the columnar defects disappears. Thus, vortex 
pinning by the columns is responsible for the upward shift 
of the superconducting transition. This is the obvious 
when one considers the shift for different columnar de- 
fect densities (matching fields B$). Figure El shows 
that a higher density of columns leads to a higher T™ ax . 

Figure|B|traces the onset temperature T^ set (H) of the 
specific heat anomaly measured in K 35 Bao.65Bi03^ii4 
along with the Tk(H)~ and Ti rr (H) lines obtained here. 
In the bismuthate compound the two latter lines coin- 
cide with the specific heat onset, which is much sharper 
than the one in YBa2Cu307-<5- This attests to the fact 
that in this compound, the effect of the columnar de- 
fects is important enough to not only remain present up 
to H C 2(T), but to actually determine the po sition of the 
upper critical field line (see also Refs. 113114]) . 

In the LLL scenario, the specific heat in magnetic field 
should behave as C p /T = F(Q), with F a universal scal- 
ing functioni22i22 Figure EJc-d) shows that the location 
in the (H, T) plane of characteristic features of C p are 
satisfactorily described by Q = constant. For example, 
vortex lattice melting in the pristine crystal occurs at 
constant Q = 7.9^ and the position of the maximum at 
Q = 3.2. However, as in Refs. I20I andl39l the magnitude 
of ACp/T = C p — Cn does not scale as expected, what- 
ever the choice for the normal contribution (Cjv) (see 
Figure Etc-d)). 



III. DATA ANALYSIS AND DISCUSSION 

A. Determination of the pinning energy 

The torque exerted by a perpendicular magnetic field 
component on the trapped vortices allows one to measure 
the energy gain due to vortex localization on the colum- 
nar defects^i i.e. the pinning energy. Writing out the 
torque step T ~ H ± M^ - H\\ M x \e=e d = -H^dG/dB^ 
shows that it is a measurement of the energy per unit 
length Ef. = $odG/dB± of vortex kinks joining different 
defects (G is the Gibbs free energy). In the limit of iso- 
lated flux lines, T w 2Ek/a\, where ao — (Qo/B) 1 / 2 is 
the flux line spacing. The kink energy Ek = (\U p ei) 1 /' 1 , 
where U p is the pinning energy of the defect per unit 
length, and £; is the vortex line tension^ For a single 
vortex line, the pinning energy can be written as^ 



U„ = 



1.2 T 2 

let 2 p 



W » 1) (3) 
(/3 «1). (4) 
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FIG. 9: YBa2Cu307-^: Closed symbols denote the irre- 
versibility lines Hi rr (T) determined from the onset of the 
third harmonic transmittivity, for single crystals with = 
(• ),B0 = 1T (■), Bj, = 2 T ( A ) and B = 5 T (V). Open 
symbols (o, □, A, V) correspond to the T™^ x -Yme tracing 
the temperature of the specific heat maximum, for the same 
crystals. The drawn lines denote the solution of Eq. If 211 . for 
matching fields = 0, 1,2, and 5 T (see text). The dotted 
line shows an evaluation of criterion choosing V = e^ b 
and a = 0.1, for B^ = 1 T only. The crossed squares depict 
the vortex lattice melting line of the pristine crystal, which is 
well described by Q = 7.9. 



The pinning strength (3 can be approximated as f3 CO re = 
(cqEo£i / kgT 2 )(l — b) 2 for vortex core pinning (for Co <C £) 
and (3 = f3 em = {^ 2 EqEi / k 2 B T 2 )(\ — b) for electromagnetic 
pinning (cq ^> The exponential factor in Eq. 0} 
expresses the reduction of the pinning energy due to vor- 
tex line wandering when /3 -C 1. As a consequence, the 
torque jump 



r = X/24-/3 1 / 2 09 » 1) (5) 

To = V^^expf-Tr//? 1 ^) (/3 «i) (6) 



is a direct measure of the pinning strength /?. Note 
that in the single vortex limit , the dimcnsionlcss torque 
plotted in Fig. Eta,b) would directly correspond to the 
square-root of the pinning strength /?. 

In practice though, experiments are rarely carried out 
in the isolated vortex limit. At high fields, only a fraction 
N t = 1 — cxp(— a^nd) ~ B^/B of the vortices are trapped 
by the columns. Then we can write 



V2 



<£o£ P 



taking (3 = l3 corei we can rewrite this as the "mean-field 
expression" 



r n = 



CO 



B d 



k B TB^' V2 2^(0)B c2 (0) 



101 



3/2 



((3 » 1) 



(7) 



(8) 

In other words, in high fields the torque jump normalized 
by the energy density scale k B T '/ 'ag£(0) is expected to 
follow an expression that only depends on the LLL scal- 
ing parameter Q. Note that the functional dependence 
on temperature and field following the LLL parameter is 
only followed provided vortex line wandering is not im- 
portant ((3 1). The exponential factor in Eq. © re- 
lated to vortex line wandering would spoil the scaling, 
because its argument cannot be expressed as a function 
of Q. The irrelevance of line wandering is slightly surpris- 
ing, for our experiments are carried out at high temper- 
atures at which it is expected to be relevant i 2 ' 4 ' 6 How- 
ever, experiment unambiguously shows that this is not 
the case. In other words, a description of pinning in terms 
of vortex line fluctuations only (i.e. the London model) 
is not a good starting point for the description of vortex 
physics in heavy-ion irradiated high temperature super- 
conductors in Tcsla fields, regardless of their anisotropy. 
The reason is that vortex line fluctuations can be inter- 
preted as the result of superposing thermally generated 
vortex loops on regular, field-generated vortices^ At 
high fields at which the LLL condition 4B£ 2 (0)/$o > Gi 
is satisfied, thermal vortices cannot be excited, as these 
imply Landau-level degeneracy^ Rather, our data show 
that the fluctuation of the order parameter amplitude 
must be considered when describing pinning in the vor- 
tex liquid, as was already suggested in Ref. |4(J 

The importance of order parameter amplitude fluctu- 
ations to pinning becomes apparent when we compare 
expression © to the experimental data, see Fig. |S{a,b). 
For YBa2Cu307_a, the experimental torque lies well be- 
low the prediction JHJ. The experimental data are consis- 
tent with either an exponential, (ao£(0) / ksT^o cx e - ^, 
or with a power-law drop, (aQ^(0)/fcsT)Fo oc Q 4 , or 
with a polynomial in Q, but shows no sign of diver- 
gent behavior on approaching the irreversibility line4£ 
In Ko.35Bao.65Bi03 however, the experimental pinning 
energy lies close to the predicted value © . Fluctuations 
have little importance here due to the small Ginzburg 
number of Ko.35Bao.65Bi03. 



B. Shift of the specific heat maximum 

We finish by showing that vortex core pinning by 
the columnar defects, such as described by Eq. |SJ leads 
to a correct estimate of the specific heat shift in both 
YBa2Cu307_,5 and Ko.35Bao.65Bi03, and may thus ac- 
count for the free energy gain obtained from pinning. We 
adopt the procedure of Ref. 0, in which the field-angle 
dependence and the defect density-dependent shift of the 
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onset temperature T^ nset {H) in K^Bai-^BiOs was well 
described quantitatively . To estimate Tg™ set (H) , we add 
the free energy change ntU p (B) from pinning to the free 
energy difference F n — F S (B) = ^u H 2 2 {^ — b) 2 of the nor- 
mal and superconducting states of the pristine material, 
and equat o 14 ! 44 



\F n - F S (B) + n t U p {B)\ = a 



k B T 
V ' 



(9) 



Here nt = N t a$ is the areal density of vortices trapped 
on a columnar defect, U p is the average pinning energy 
per vortex per unit length, given by Eq. i jffi . V = e£ 3 
is the coherence volume, and a < 1. In Ref. Il4l 
Eq. (PJ was found to quantitatively reproduce all lines 
\T irr (H), T k {H),T^ set (H)} if one chooses U p to be given 
by the expression for electromagnetic pinning, Eq. (@J) 
with (3 = P em and awl. There is some latitude in 
the choice of parameters: the dotted line in Figure [B] 
shows that similar good agreement can be obtained with 
f3 = f3 core and a = 0.1. 

We find that the position of the T™ ax -Yme in 
YBa 2 Cu 3 07_,5 can also be described in this way. The 
dotted line in Fig. 01 shows the solution of Eq. OH with pa- 
rameters for YBa2Cu307_,5, the pinning energy as given 
by Eq. (@J with f3 = /3 core , and a = 0.1. 

Another approach is to choose the correlation volume 
V = £<2q£ rather than e£ 3 . In the absence of columnar 
defects (n t = 0), Eq. © can then be reduced to 



Q 



3/2 



47rv / 2a, 



(10) 



which, for a = 1.25 (Q = 7.9), perfectly describes the 
position of the vortex lattice melting line (see Fig. [§)). 
In the presence of columnar defects, one again adds the 
free energy gain n t U p . Taking the same expressions U p = 
(c /2£) 2 £ (l-fr) 2 and n t = -B^/<£>o used to derive Eq. ©, 
Eq. (jHJ) becomes 



)3/2 



1 + B^/S c2 (0) 



(11) 



In other words, the criterion © takes the form "Q = 
constant" , where the constant depends on the density of 
columnar defects. In order to describe the shift of the 
specific heat maximum with B^, we adopt the critical 
temperature T^ IF = 93.1 K obtained from the scaling of 
the torque and the specific heat, and evaluate 



Q = 3.2/[l + 5 /£ c2 (O)] 



2/3 



(12) 



so as to recover the correct position of the specific heat 
maximum for B$ = (see Fig.|HJ). The dependence l|12|) 
well describes the upwards shift of the specific heat max- 
imum in a magnetic field (Fig. |5J). Thus, we find that 
the incorporation of the pinning energy in the free en- 
ergy difference between the normal and superconducting 
states satisfactorily describes the evolution of the specific 
heat as function of defect density. 



IV. SUMMARY AND CONCLUSIONS 



Thermodynamic measurements in the vortex liquid 
of heavy-ion irradiated YBa 2 Cu 3 07_,5 show that, as in 
K x Bai_a;Bi03, the reduction of the average free energy 
of the superconductor due to vortex pinning on columnar 
defects remains important all the way into the regime of 
strong order parameter amplitude fluctuations. Notably, 
the effect of pinning on the free energy is sufficient to 
shift the superconducting transition, as measured by the 
heat capacity, further upwards as the defect density in- 
creases. The scaling of the pinning energy with the LLL 
parameter Q shows that pinning is affected by fluctua- 
tions of the overall order parameter amplitude. These 
should therefore be taken into account in any descrip- 
tion of pinning in the vortex liquid - a model based on 
vortex line positional fluctuations only is bound to be in- 
adequate. However, due to the lower Ginzburg number, 
order parameter fluctuations are much less important in 
the bismuthate superconductor than in YBa 2 Cu307_,5. 

This brings us to the final point, which is the Bosc- 
glass transition. Our results indicate that depinning of 
vortices from columnar defects due to line wandering is 
not likely to be a very effective mechanism for this transi- 
tion. On the other hand, fluctuations of the overall order 
parameter amplitude are also expected to become unim- 
portant at, or near to, the Bose-glass transition. We 
note that, contrary to the melting line in the pristine 
crystals, the irreversibility line (and especially its high- 
field part) cannot be described as a line of constant Q 
and therefore does not satisfy LLL scaling. We therefore 
suggest that another derealization mechanism drives the 
Bose-glass transition. This could be, for example, vortex 
derealization due to variations of the pinning potential. 
Analogously to the electronic mobility in disordered semi- 
conductors, the most weakly bound vortices could be de- 
localized at a much lower temperature than the bulk of 
the vortex matter £ Another possibility is that the plastic 
properties of the vortex ensemble play a much more im- 
portant role than hitherto considered. For example, the 
trapped vortex ensemble could constitute a polycristal, 
vortex derealization taking place on the grain bound- 
aries. In both cases, another lower energy scale than 
the average pinning energy is involved in vortex dereal- 
ization at Ti rri a scale that manifests itself through the 
large separation between T irr and [Tfc(iJ), Tc p (H)] lines 
in YBa2Cu 3 07_5. 

The situation in cubic K^Bai-^BiOa seems to be dif- 
ferent. We have shown that the disappearance of pinning 
by columnar defects in this compound happens indistin- 
guishably close to the superconducting-to-normal bound- 
ary. It means that in this compound, vortices are neither 
depinned, nor delocalizcd at all, and the Bosc glass phase 
subsists up to the temperature at which superconductiv- 
ity disappears altogether. 
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